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Electromagnetic polarizabilities of the nucleon are analyzed in a hedgehog model 
with quark and meson degrees of freedom. Semiclassical methods are used (linear 
response theory, quantization via cranking). It is found that in hedgehog models 
(Skyrmion, chiral quark models, Nambu-Jona-Lasinio model), the average electric 
polarizability of the nucleon, a at, is of the order N c , and the splitting of the neutron 
and proton electric(proper) polarizabilities, 5a = a n — a p , is of the order 1/N C . We 
present a general argument why one expects 5a > in models with a pionic cloud. 
Our model prediction for the sign and magnitude of 5a is in agreement with recent 
measurements. The obtained value for a^, however, is roughly a factor of three too 
large. This is because of two problems with our particular model: a too strong pion 
tail, and the degeneracy of N and A states in the large- N c limit. This degeneracy 
also results in a very strong ^-dependence of the paramagnetic part of the magnetic 
polarizability, (3, which is of the order N%. We compare the large- results to the 
one-loop chiral perturbation theory predictions, and show the importance of A-effects 
in pionic loops. We also investigate the role of non-minimal substitution terms in the 
effective lagrangian on the polarizabilities of the nucleon. 
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I. INTRODUCTION 



Polarizabilities are important fundamental properties of particles - they determine dy- 
namical response of a bound system to external perturbations, and provide valuable insight 
into internal strong-interaction structure. Recent measurements of electromagnetic polariz- 
abilities of the proton and neutron H renewed theoretical interest in these quantities 
||. Many calculations of the electric, a, and magnetic, (3, polarizabilities of the nucleon 
can be found in the literature in models ranging from various quark models and bags |§-|3| 
cloudy bags |T(| to Skyrmions [PH13|. A quenched lattice calculation for a has also been 
done |T4[]. Extensive reviews of the subject are available ]To]JTB| . 

In this paper we analyze the nucleon polarizabilities in hedgehog models, using the 
framework of the linear response theory. Our analysis differs significantly from previous 
works by allowing the soliton to deform ( dispersive effects). Our semiclassical methods 
are described in detail in the following paper ||17|| , referred to as (II). The reader who is not 
familiar with the basics of hedgehog models or semiclassical methods used in their treatment, 
is urged to read (II) before this paper. 

In recent years numerous hedgehog models, such as the Skyrme model [IB1-BI], chiral 
quark models p2|-p7|1 , hybrid bag models [p8] , chiral models with confinement [ p6| , p9| -|31| , or 

P5|-PD|], were quite success- 



the Nambu-Jona-Lasinio model in the solitonic treatment 
ful in describing the phenomenology of nucleon structure. The basis for these hedgehogs 
models is the large-iV c limit |40| , |4lH of QCD. Thus, l/N c should consistently be treated as 
an expansion parameter in this approach. Previous analyses of polarizabilities in hedgehog 
models 



11-13 



have not emphasized the need for consistent N c - counting. In particular, a 
frozen approximation in the treatment of the response to electromagnetic perturbations was 
made, and the effects of distortions of the soliton, which occur at a relevant iV c -level, were 
not included. In this paper we develop consistent iV c -counting rules for a and (3, and eval- 
uate these observables in a specific model. We point out serious difficulties arising in the 
case of the magnetic polarizabilities. Special attention is paid to the issue of neutron-proton 
splitting of the polarizabilities, which has already been reported elsewhere |^2| . 

Throughout this paper we use the so-called proper (or static) polarizabilities, a and (3, 
defined via the energy shift of an object in constant electric, E, or magnetic, H, fields |[43|| : 



5£ 



-\*E* 



These are directly related (see Sec. |l I FF] ) to the coefficients a and (3 encountered in the 
Compton amplitude plH47|Jl5|JT6|JT^] 



a 



a + Aa, Aa 



3M 



/3 = /3, 



[1.2) 



where the recoil term Aa involves the charge, Q, mass M, and electric mean square radius 



of the particle. Although the experimental errors are still substantial P5| , the new measure- 
ments indicate that the proper electric polarizability is larger for the neutron than for the 
proton, 



Oir, 



(12.0 ± 1.5 ±2.0) x 10~ 4 fm 3 , 
(7.2 ± 1.0 ±1.0) x 10~ 4 fm 3 , 



;i.3) 
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and that the magnetic polarizabilities are positive, and equal within experimental errors: 



/3 n = (3.1 =p 1.5 =f 2.0) x 10~ 4 fm 3 , 

p p = (3.4=f1.0=f1.0) x fm 3 . (1.4) 

A number of results presented in this paper are generic to any hedgehog model, e.g. the 
iV c -counting rules and the role of the intermediate it — A states. Our numerical results are 
obtained using the simple model of Ref. |22| (Sec. II of (II)). Since this model possesses 
both quark and meson degrees of freedom, it is in this respect representative of a general 
class of hedgehog models. Our calculations can be repeated in any hedgehog model with 
minor modifications. 

The organization of this article is as follows: In Sec. || we write down the gauged er- 
model lagrangian, introduce collective coordinates in the usual way (Sec. Ill F of (II)), 
and classify various electromagnetic perturbations according to the hedgehog symmetries. 
We discuss the sea-gull and dispersive contributions to the Compton amplitude (Sec. |II B| ) . 



and present the iV c -counting rules (Sec. [11 Dp. For the electric polarizabilities we find that 



the neutron-proton splitting effect is two powers of N c suppressed compared to the average 
nucleon value: 

«jv = («„ + ot p )/2 ~ N c , 

5a = a n — a p ~ 1/N C . (1.5) 
For the magnetic case the ^-counting is more complicated (Sec. |TV| ). 



In Sec. 1 1 J we obtain numerical results for the electric polarizabilities. For a at, the sea- 
gull contribution (Sec. | LI A| ) is dominant, and the valence quark effects enter at the level 
of 10 %. This is due to the long-range nature of the pion. Our model prediction is about a 
factor of three too large than the experimental number. We discuss this discrepancy, which 
is due to two resons. Firstly, our specific model has a pionic tail which is too strong (i.e. 
the value of g n NN is too large compared to nature). This enhances the model prediction 
by about a factor of 2. A more fundamental reason is discussed in Sec. [V], where we 
point out that hedgehog predictions for some observables largely overestimate results due to 
the implicit treatment of the A resonanse as degenerate with the nucleon. Our prediction 
for 5a = a n — a p are less sensitive on the strength of the pionic tail, and they are not 
affected by problems discussed in Sec. |V|. Numerically, we obtain 5a = 5.4 x 10~ 4 /m 3 [J2 



in agreement with the sign and magnitude given by the experiment. This splitting arises 
naturally in hedgehog models when dispersive effects are taken into account. The effect is 
dominated by the distortion of the pionic cloud. We also present a classical argument why 
the sign of 5a is expected to be positive in models with pionic clouds (Sec. |IIE|). 



In Sec firT], we show that the standard retardation correction (Eq. |1.2| ) strictly holds 



in our linear response method. In Sec. [IV| we discuss the issue of magnetic polarizability in 
hedgehog models. We show that the degeneracy of the A and nucleon masses in the large- N c 
limit precludes the use of linear response theory to determine (3. The N — A paramagnetic 
contribution, /3tva (given by the Born term with the intermediate A state) dominates the 
iV c behavior of (3, and leads to 

Pn = (f3 n + (3 p )/2~N 3 , (1.6) 
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In our view, this invalidates the claims of Refs. [pTT^ [T3[| that the diamagnetic sea-gull inter- 
action cancels the paramagnetic A term, since these terms occur at different iV c -levels In 
Sec. [IV C] we also show the relevance of dispersive contributions in the calculation of the 
magnetic polarizability. 

In Sec. |V| we compare our results to predictions of the chiral perturbation theory (%PT). 
We find, that in the chiral limit the hedgehog predictions and the xPT predictions for a 
and (3 agree up to a factor which can be attributed to the role of the A in pionic loops. 
We discuss interesting physical implications of this issue, and point out that while the naive 
hedgehog predictions overestimate the role of the A resonance in hadronic loops, the naive 
approach to xPT drops these important contributions altogether. We discuss the effects of 
finite N — A mass splitting in these approaches, and propose how the hedgehog and xPT 
results should be corrected. 

In Sec. |VI] we study the effects of the pionic substructure on nucleon polarizabilities. 
These effects are introduced via non-minimal substitution terms £9 and C w [^,0] in the 
effective lagrangian. At the mean-field level, we obtain simple expressions involving pion 
polarizability and pion mean squared radius. Numerically, these pion structure contributions 
enter at the level of ~ 1 — 2 x 10~ 4 /m 3 , which is a small but non-negligible effect. 

In Sec. |V11| we make a few comments relevant to calculations in other hedgehog models. 



Throughout the paper we use the units h — 1, c = 1, e /(he) = 1/137. 



II. ELECTROMAGNETIC PERTURBATIONS 

Our basic methods of treatment of external perturbations in hedgehog models are de- 
scribed in (II). The approach developed in (II) is a rather straightforward adaptation of the 
RPA method of traditional many-body physics. The diference with the traditional nuclear 
physics case is the presence of mesonic degrees of freedom, relativistic dynamics, and special 
symmetries (Sec. II B in (II)). 

Polarizabilities are defined as second-order energy shifts due to (static) external elec- 
tromagnetic fields, Eq. ( |1.1|) , or, equivalently, via Compton amplitude (Sec. |II B|) . They 



arise from two types of terms: sea-gulls, which are due to quadratic terms in E or B in the 
lagrangian, and from dispersive terms. Sea-gull effects result in local expressions involving 
the pion mean field [|TT|-|T3"|, while the dispersive terms are given by the usual second-order 
perturbation theory expressions, which result from linear perturbations in the lagrangian. 
These dispersive effects distort the soliton. They are particularly important in the neutron- 
proton splitting effects, as well as in the magnetic response, and have not been considered 
in earlier works. More importantly, these effects enter at the same A c -level as the sea-gull 
terms, and thus should be included in order to comply with the basic organizational principle 
of the hedgehog approach, namely, the l/A c -perturbation expansion. The dispersive terms 
in the Compton amplitude correspond, in our approach, to diagrams in which the propa- 
gator between the two photon insertions is the RPA (or linear response) propagator (Fig. 

[I] ). Relevant expressions are obtained using equations-of-motion method, which amounts 
to solving linear differential equations with potentials and driving terms determined by the 
solitonic profile. In addition, electromagnetic perturbations are carried in the presence of 
cranking, which ensures projection on states with good quantum numbers (Sec. II G in 

(II) ). 
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A. Lagrangian in presence of electromagnetic interactions 



The first step is to identify perturbations resulting from coupling to electromagnetic 
interactions, which are introduced by gauging the cr-model lagrangian (II.2.1): 

C = $ [i@ + g (a + Z7 5 r • 7r)] if? 

+ \{d»af + |(^7t) 2 - U (a, tt) . (2.1) 

Through this procedure one generates a dispersive interaction, as well as a covariantizing 
sea-gull term: 

C disp . = ~eA^ fa (-^- + §r 3 ) 7^ + x ^ 7r ) 3 ) > 

£ S ea- g u« = |e 2 A^ (TT 2 + 7T 2 ) , (2.2) 

where ^ is the quark field, 7r is the pion field, and A^ is the photon field. It is hoped that 
this minimal substitution procedure gives the bulk of electromagnetic interactions in our 
effective hadronic lagrangian, however, non-minimal substitution terms may a priori play an 



important role. In Section [V| we discuss the role of simple non-minimal substitution terms 
in effective lagrangians, and find that their effects enter (in the electric polarizability) at the 
level of 10 %, compared to the terms resulting from (|2.2| ). 

To identify the appropriate perturbations, we make a transformation of the full la- 
grangian (II. 2.1, |2.2j ) to the isorotating frame, according to Eqs. (II. 3. 22). This is done 
because we are interested in the linear response of the nucleon, which is "projected out" of 
the hedgehog (see Sec. Ill in (II) for details of projection via cranking, and linear response 
in presence of cranking). For a more compact notation, we also replace the quark field 
bilinears by N c times the bilinears of the valence quark, q. The Lorentz gauge is used, and 
A = -r ■ E, A = -\r x B. We obtain 

C — > C° + C\ + C\ + C\o + C\i + £jgi + £ 2 £i 

+ ~\- Cj^i -\- £^i, (2-3) 

C° = N c q(i$ + g(a + ij 5 T ■ tt)) q 

+ \{d^af + i(^7r) 2 - U (a, tt) , (2.4) 

4 = -iWQA-T)g-A-(7T xtt), (2.5) 
C\ = \{\ x tt) 2 , (2.6) 

C\ = er-E^q, (2.7) 

C x E i = er ■ E (N e q^c - rq + c - (ir x tt)) , (2.8) 

£% = ie 2 (r-£) 2 (cx7rf, (2.9) 

C\ E1 = - er ■ E(c x tt) • (A x tt), (2.10) 
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4*o = -ie(rxB)Va<g, (2.11) 
^si = -|e (r x B)*(A^ C q ] \c ■ ra { q 

- c-(tt x Viir)), (2.12) 

£|i = -^e 2 (r x B) 2 (cx tt) 2 , (2.13) 



where A is the cranking velocity, and the collective vector c = |Tr[r 3 i?T-BT] is defined in 
App. B in (II). In the above equations, subscript A refers to term arising upon cranking, E° 
and E 1 denote isoscalar and isovector electric perturbations, and B° and B 1 denote isoscalar 
and isovector magnetic perturbations. Indices i in Eqs. (|2.11| - |2.12|) are spatial Cartesian 



indices. Superscript 1 denotes dispersive terms, which lead to linear shifts in fields according 
to perturbation theory. Superscript 2 denotes terms which are quadratic in perturbations. 
They include a term quadratic in A (Eq. [2.6|) , sea-gull terms (Eqs. p^9|J2.13|) , and a mixed 
term (Eq. |2.1(J ), which will play an essential role in the splitting of the neutron and proton 



polarizabilities. The classification of various dispersive terms, as well as the explicit forms 
of the corresponding sources arising in linear response equations, are given in Table |. The 
grand spin (sum of spin and isospin) is denoted by K, parity is denoted by V, and the 
grand-reversal symmetry by TZ (see Sec. II B in (II) for details). 



B. Compton amplitude and polarizabilities 

The Compton amplitude corresponding to lagrangian (|2.3j ) reflects the presence of both 
dispersive and sea-gull terms. It can be written as [12] 

M„ v (p, q) = M^(p, q) + S„ v (q), 
M^ip,q) = 

i J <fixe**W)\T (j;- m ix)r u m iO)) \N(p)), 

(2.14) 

where M M „ is the dispersive T-product part, and S^ u is the sea-gull part. Correspondingly, 
the polarizabilities have dispersive as well as sea-gull parts, and can be written as fll2"| 

& &disp. & sea— gull} ft fidisp. Psea— guilt 

\(N\Jd 3 xE-rJ*- m ir)\b')\ 2 



2E 



Ey — / - \ 



Oisea-gull = {N\ J d 3 X {E ■ r) 2 S 00 \N), 

\(N\$d 3 x\B-(r x J e - rn -(r))\b')\ 2 



ftdisp. 2 ^ ' 



y Ey — En 



Psea-guii = \(N\ J d 3 x e mni B m r n e klj B k r l Sij\N) . 

(2.15) 

Hats denote unit vectors in the direction of E or B fields, and \b') is an intermediate state 
with energy Ey. 
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The sea-gull parts of expressions ( |2.15|) can be readily identified from Eq. ( |2.9j , p.l3|) . 

One gets 

S 00 = e 2 (c x tt) 2 , 

S tj = -e 2 5 ij (c x tt) 2 . (2.16) 

Evaluating these in the collective nucleon state one obtains (App. B in (II), p!2|JT3| ) 

87re 2 f 

asea-guii = -2f3 sea „ guU = — — J dr r 4 ir 2 h , (2.17) 



where ir^ is the hedgehog pion field profile. In Eq. ( 2.17 ) we recognize the hedgehog 
model relation between the seagull contributions to the electric and magnetic polarizabilities 
12 , 13] . The radial integral in Eq. ( |2.17 ) is quadratic in the pion field. In the Skyrme model 



additional terms arise due to higher-order terms in the lagrangian DTTHTB 



C. Linear response 



In our treatment of dispersive contributions, the intermediate states \b') in Eq. fl2.15|) 
correspond to RPA excited states, and the energy denominators involve energies of these 
excitations (Fig. [I]). These contributions can be written in the general form given in Eqs. 
(II. 3. 27). The electric and magnetic dispersive polarizabilities result from either isoscalar 
dipole or isovector dipole transitions, and will be labeled by E° and E 1 , or B° and B 1 , 



respectively. These dispersive contributions, as well as the sea-gull contributions ( |2.17D 



contribute equal amounts to the neutron and the proton. This is because the resulting 
collective operators in Eqs. (II. 3. 27) are isoscalar (bilinear in the collective vector c), and 
matrix elements are equal for the proton and the neutron. The neutron-proton splitting 
effects arise from the mechanism described in Sec. Ill G in (II). This can be briefly described 
in the following way: We perform electromagnetic perturbations not on hedgehogs, but on 
nucleons. These are obtained from hedgehogs via projection (cranking), which in turn may 
also be viewed as linear response. It is sufficient to work to linear order in the cranking 
velocity, A, and we can write down symbolically 

\N) = (l + gV cr )\H), (2.18) 

where Q is the RPA propagator, V cr is the cranking interaction, and \H) is the hedgehog 
state. Using these "cranked" states in our perturbation theory leads to expressions of the 
form of Eqs. (II. 3. 29), with one cranking, and two electromagnetic interactions. As discussed 
in Sec. Ill G in (II), the K v numbers of interactions in Eqs. (II. 3. 29) have to be additive to 
K v = + , since the hedgehog has K v = + . From Table | we can see that we can compose 
cranking, one isoscalar electric (magnetic), and one isovector electric (magnetic) interaction 
to K v = 0+. Since V cr carries a collective operator A, the appropriate collective operator is 
isovector, and neutron-proton splitting of polarizabilities is generated. 

From a slightly different but equivalent point of view, we may understand the isospin 
effect in electric polarizabilities by inspecting the term ( |2.10| ). It is linear in cranking, and 
linear in isovector electric {E 1 ) perturbation. It leads to the following contribution to a: 
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a mes. = _ e ( N \ J d 3 r(E-r)[(cx tt) ■ (A x 7r)]|iV). (2.19) 



mes. 
A 



If the pion field in the above equation were taken to be just the hedgehog profile, then a 
would vanish by parity. However, the isoscalar electric perturbation (Section [III C ) distorts 



the meson fields, and 7r = + 57r E o, where 8tz e g has S- and D-w&ve components. As a 
result, the integral is non-zero. The corresponding collective operator is c - A = I3/O, where 
is the moment of inertia (see App. B in (II)), and we obtain 

a rnes. = _ g h g/ g [ @ . r )(lT h ■ 5l7 E o), . (2.20) 





From the quark parts of Eqs. (|2.3|) we get 



+ 5q E0 \eE ■ rc ■ T5q x }\N) + h.c, (2.21) 

where various terms in the interaction are sandwiched with shifts corresponding to other 
interactions. Eqs. (|2.20| , |2~2l|) can be straightforwardly obtained from Eq. (II. 3. 29). Terms 
( j2.20j , [2.2lD are proportional to the nucleon isospin, J 3 , thus are responsible for splitting of 



the neutron and proton electric polarizabilities, 5a. 

For the splitting of magnetic polarizabilities, mesons do not contribute at the linear- 
response level, since there is no mesonic analog of the term ( |2.10|) in lagrangian ( |2.3j ). Only 
the magnetic analog of the quark part (Eq. |2.21| ) is present. However, the linear response 
calculation of the magnetic polarizability encounters fundamental problems, which will be 
discussed in Sec. [TV]. 



D. N c - counting 

In Eqs. |2.3| we list explicitly occurrences of iV c . We also recall ETJ that in the large- N t 



c 



limit the meson fields scale as yN cl the quark fields scale as 1, and the moment of inertia 
scales as N c . Note that the sources for isoscalar electromagnetic perturbations are one power 
down compared to the corresponding sources for the isovector interactions. We easily arrive 
at the following iV c -rules for electric polarizabilities: 

otsea-guii ~ N c , a E o ~ l/N c , a E i ~ N c , a x ~ 1/N C . (2.22) 
This leads directly to the following rules for the nucleon polarizabilities: 

(X n + CYp 

= ~ = a se a-gull + Ot E l ~ N c , 

8a = a n — a p = a\ ~ 1/N C . (2.23) 

Note, that since the a E o part of the dispersive term contributes to ajy at a subleading 
level, it should be dropped according to A^-rules. Other physical effects affect our results 
at this level, e.g. centrifugal stretching, center-of-mass corrections to the soliton mass, 
^-suppressed terms in the effective lagrangian, etc. Numerically, the a E o contribution is 
negligible, confiring the validity of ^-counting for polarizabilities. 

Naively, one would write down expressions analogous to Eqs. (|2.23|) for the magnetic 



polarizability, (3. We will show in Sec. |V]that this is not correct. 
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III. ELECTRIC POL ARIZ ABILITY 



In this section we present our numerical results for a. Numerical methods, explicit forms 
of the differential equations, and other details are given in (II). The model parameters for 



our numerical calculations are taken from Refs. 



93MeV, gF n = 500MeV, 



139.6MeV, and m a = 1200MeF. 



A. Sea-gull contribution 

Because the pion field is long ranged, the nucleon electric polarizability, a at, is dominated 
by the pionic sea-gull contribution (|2.17| ). Numerically, in our model we get 



cwfluii = 28 x 1(T 4 fm 6 , (3.1) 

which is roughly three times larger than the experimental value in ( |1.3|) . There are two 
basic reasons for this discrepancy. The first one is discussed below, and the other in Sec. 
|V[ Figure (0) shows the radial density of the integrand of ( [2.17| ) (solid line). It is clear 



that the sea-gull term acquires most of its value from the asymptotic region, r > lfm. 
Recall, that the asymptotic behavior of the pion tail in hedgehog models is determined by 



the pion-nucleon coupling constant, g^NN P0|j51| - Using the Goldberger-Treiman relation, 
we can write down 



7T 



asympt. _ ^ g \%^F v )x{m^ + 1/ 'r)exp(-m^r) / V. (3.2) 



Our model predicts g& = 1-86 [[SlU , hence the tail contribution is overestimated by a factor 
of (gT del / 9a P T ~ 2 - The dotted line in Fig. (@) shows the integrand of Eq. (|2~T7| ) with the 
pion field having the form (|3.2j ) in the whole region of r. This corresponds to the chiral limit 
case, and will be discussed in detail in Sec. M. We can see from Fig. (fj) that the size of 
a sea- g uii is controlled by how fast the solid line departs from the dashed line at lower values of 
r. In other words, the results are sensitive to the profile of the pion field in the intermediate 
region of about lfm. Note, that 7T/j enters the expression (|2.17|) quadratically. Therefore, 
we expect substantial model sensitivity in hedgehog predictions of a S ea-guih e.g. models 



which have a suppressed pion field, such as hedgehog models with confinement [[26| ,p9|-|3T 
or models with vector mesons, e.g. |52"1J53~1| , are expected to predict lower values for a sea - gu ii- 
We stress that to compare fairly the model predictions with experiment, the model should 
predict correctly the quantity gA/F n , which enters the asymptotic form (|3.2|) . 

In Sec. [V] we show how the sea-gull contribution is additionally suppressed when N — A 
mass splitting effect is taken into account. This results in another factor of ~ 2 reduction. 
With these corrections we note that it is possible to put the sea-gull contribution in the 
right experimental range. The model uncertainty, however, is big. 



B. Isovector electric perturbation 

In Sec. lirOwe have shown that the dispersive electric isovector (E 1 ) contribution to 
polarizability, a^i, is of the order iV c , the same order as for a sea ~guii- The E 1 perturbation 
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in lagrangian ( |2.8|) has K vn = , 1 and 2 (Table |). Since it is odd under grand- 
reversal 1Z, only the valence quark fields acquire shifts (Sec. Ill A in (II)), and we solve 
equations of the form of Eqs. (II. 3. 7): 

(h- e)q E i = \N c er ■ Ec rq h . (3.3) 

These equations are decomposed into K = 0, 1 and 2 components, as described in App. A 
in (II), and solved numerically. The corresponding polarizability is calculated according to 
the expression 

a E i = —2N c e (N\ J <?xq E r ■ Ec ■ rq h \N), (3.4) 

where \N) is the nucleon collective state (App. B in (II)). Note that with our definition, the 
spinor q E i carries the collective variable c, such that Eq. ( |3.4|) contains a matrix element of 
a collective operator quadratic in c, similarly to the case of the sea-gull part in Sec. |II B| . 
Numerically, we get in our model 

a E i = 3.49 x 10~ 4 fm 3 , (3.5) 

with K = 0,1, and 2 components of q El contributing —12%, 26%, and 86%, respectively. 
The dashed line in Fig. ^| depicts the radial density of a E i. Because the quark mass in 
our model is 500MeV, the purely quark a E i contribution is strongly suppressed compared 
to the pion sea-gull contribution. However, the value of a E \ is non-negligible compared to 
experimental numbers (Table |I]), and quark effects are substantial. 

In models with vector mesons |)2]|k|, the E 1 perturbation involves shifts in the p and uj 
mesons, and these effects must be included in calculations of the electric polarizability ||54| . 



They enter at the leading-iV c level, and should contribute comparably to the valence quark 
effects. 



C. Isoscalar electric perturbation 

Now we turn to a more difficult case of the isoscalar electric perturbation, E°, which 
has K vn = 1~ + . These are quantum numbers of the translational zero mode, and, as 
explained in Sec. Ill A in (II), we have to deal with full RPA equations of the form (II. 3. 3). 
The appearance of the zero mode is easy to understand. The hedgehog soliton possesses 
isoscalar electric charge, Q I=0 = \e. Thus, in a constant electric field, the soliton accelerates 
in the direction of E — the translational mode is excited. In addition to this zero-mode 
motion, the soliton is getting distorted (as viewed from the center of mass). Below we 
show that the excitation of the translational motion corresponds to the Thompson limit of 
the Compton scattering, and the deformation corresponds to the polarizability. In order 
to separate the zero mode from the physical modes, instead of a constant electric field we 
consider a slowly time- dependent field of the form 

E{t) = \ (Eve** + E*e iujt ) (3.6) 

The resulting sources in Eq. (II. 3. 3) are 
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Jx=Jy = \eN c 1 (r-E)q h) 

Jz = Jp = 0. (3.7) 

Using the technique described in Sec. Ill B in (II), we find that the total isoscalar electric 
polarizability consist of the translational zero-mode part, and a "physical" part: 

„tot. „zero i ^phys. fo q\ 

a E o = a E o + a E o , (3.8J 



where the zero-mode part diverges in the limit ui — > 



(Q 



I=0\2 



and the physical part, a p E o S \ has a finite u — > limit. The quantity M in Eq. ( P-9|) 
is the soliton mass. In the expression for the forward Compton scattering amplitude, a 



is multiplied by u 2 H47lJ44]J16| . We note immediately that uj 2 a E S° = — ^ M - is just the 



Thompson term in scattering of a particle with charge Q 1 and mass M. Thus, the zero- 
mode part is responsible for the Thompson limit. Since RPA leads to small fluctuation 



equations of motion, analogously to the case of classical physics ||44J1 , it is clear it leads to 
the correct Thompson limit. 

In the case of the E l perturbation, the quantum numbers precluded excitation of a zero- 
mode on top of the hedgehog solution. This is a manifestation of the fact that the hedgehog 
does not have any isovector electric charge, Q I=1 = 0. This charge arises only upon cranking, 
but the organization of our perturbation theory is such, that we treat cranking, E° and E l 
perturbations separately (Sec. Ill G and V in (II)). If linear response were performed on a 
cranked soliton, then the full charge Q = Q I=0 + Q I=1 would appear in Thompson scattering. 



Earlier works ||12| , |T3l| on electric polarizabilities in hedgehog models did not take into 
account the effects of zero modes. In fact, this is justified by the 1/N C expansion. The E° 
contribution is suppressed by two powers of N c compared to the E 1 contribution (Eq. ( [2.22 ). 



Therefore, if one is interested in the leading N c behavior of a at, it is sufficient to consider 
the isovector electric perturbation only, where the issue of the translational zero mode does 
not arise. 

Since we are interested in the splitting of the neutron and proton polarizabilities, we 
have to calculate the quantities a™ es - (Eq. ( |2.20| )) and a q ^' ark (Eq. ( |2.21| )). Thus, we have 
to find the shifts in the fields due to the isoscalar perturbation, E°. We need to extract the 
"physical" parts of the solution. The numerical procedure has been described in Sec. Ill 
C in (II). Here we only remark, that very good numerical accuracy is necessary in order to 
separate the zero mode from the physical mode. This is because at small values of u> the 
solution can be written as 

iathQ Q 

^ = " ~Mu + W + (3 - 10) 

where £o is the translational mode, £i is the conjugated mode ("boost mode"), and £ p h ys . 
is the physical mode which we want to extract (Sec. Ill B in (II)). Solutions of equations 
(II. 3. 3) in the limit uj —>■ give full £, from which we subtract the pieces divergent in u. 
This can easily be done, since we know the exact forms of £o and £i (App. C in (II)). The 
presence of the zero mode provides a useful algebraic check of equations from App. A in 
(II), since the divergent parts of Eq. (|3.10| ) have known coefficients. 
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D. Neutron-proton splitting of electric polarizabilities 



Having solved the linear response equations for the E° perturbation, we use q p ^o S ' = 
XP h v s -+YP h y s and phys - = 2Z ph y s - (App. A in (II)) in equations (gg(| ggj] ). The dominant 
part to a\ comes from the pionic contribution. The pion shift has S- and D-w&ve components 
(Table |), which we denote by tts and no, respectively. The explicit expression for the 
mesonic contribution of 5a mes has the form 




The quark contribution is obtained analogously from Eq. ( [2.2 1| ). The numerical results give 

5a = -5.6 x 1(T 4 fm 3 , (3.12) 

with the quarks carrying —5% of the total. Hence, as in the E 1 case, we observe the 
dominance of the pionic contribution. Figure [3] shows the radial density of expression ( ^.11| ). 
The long-range nature of the pion is evident. As discussed in Sec. [HI B| , our predictions are 
enhanced due to the fact that g& has a too large value in our model. In the present case, 
however, we estimate this effect at the level of g™ odel j g^ p ' ~ 1.4, rather than a factor of 2 in 
the pionic sea-gull term. This is because Eq. ( |3.11|) is linear in Kh (recall that asymptotically 
iih is proportional to Eq. (|3~2|)), and its and no are not dependent on the strength of 
the pionic tail in TTh- The reason is that the sources in the linear response equations come 
from the quarks, which are short-ranged. In the asymptotic region, where Eq. ( |3 . 1 1| ) gets 



most of its contribution, the shifts in the pion field, tt$ and Tin, depend on the total strength 
of the source (Q I=0 — |), and the pion mass, which enters the Green's function, but not on 
the strength of the tail of 7T/j. Consequently, our overestimate of 5a should be at the level 
of 40% only. 

Our results are summarized in Table [H[ We note, that the sign of the neutron-proton 
splitting effect, and its magnitude, agree with the recent experimental numbers. 

E. Pionic cloud and sign of 5a 

In the previous section the sign of 5a = a n — a p was found to be negative. The following 
plausibility argument can be given as to why we expect this behavior in models with pionic 
clouds. The hedgehog models imply that the nucleon consists of a quark core, carrying 
isoscalar charge, and a pion cloud, carrying isovector charge. In the proton, these charges 
have the same sign, and the electric field distorts both the core and the cloud, but it does not 
displace them relative to each other. In the case of the neutron, the cloud and the core, in 
addition to being deformed also get displaced, since their charges are opposite. This results 
in additional polarizability. Calculations based on other methods also give a n > a p |53| , |5*U 



F. Retardation term 

Our calculation was performed with a spatially constant electric field. Now consider 
plane-wave photons, as in Compton scattering, i.e. 
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E = E exp(ik± ■ x) 



E (l + ik± 
u 2 . 



X 



l(fcj 



x 



(3.13) 



The linear-response source, j, is modified accordingly, and the charge Q = (Sec. Ill in 
(II)) is replaced by an effective charge Q e ^' = Q(l — \uj 2 {t 2 ) E ). Substituting this expression 



to the expression for the zero-mode part of the polarizability, Eq. ( |3.9[ ), we obtain the usual 
retardation term, exactly as in the classical derivation of ref. ||44|| , as well as in more general 
derivations [ffBlHTl. 



IV. MAGNETIC POLARIZABILITY 

Although naively one would think that the analysis of the magnetic polarizabilities would 
parallel the electric case, it turns out that there are fundamental difficulties. The reason 
is the non-commutativity of the u — > limit of the Compton scattering, and the large- N c 
limit, in the magnetic case. As stressed earlier, since the large- N c limit is the basic principle 
behind the hedgehog approach, one has to work consistently in the N c counting. We point 
out that with the linear response approach to the magnetic polarizabilities it is not possible. 
We also show importance of pionic dispersive terms in the magnetic polarizability, which 
enter at the same level as the pionic sea-gull contribution. Such terms were neglected in 



previous works [ITTl-O] 



A. The N - A Born term 

The importance of the Born term with intermediate A state in the estimates for (3 is a 
well recognized fact. The hedgehog possesses isovector magnetic moment, and the magnetic 
interaction term in the effective hamiltonian has the form flTT| , |T3"| 

E magn - = 3// =1 c ■ B, (4.1) 

where // =1 is the isovector magnetic moment, and the collective vector c is defined in App. 
B in (II). Hamiltonian ( |4.1|) leads to iV — A transitions, and, according to perturbation 
theory, the contribution to the magnetic polarizability is 



M A - M n M A - M 



N 



where in the last equality we have used Eq. (II. B. 9). Now, recalling the Ac-counting rules 
T9| , [5ll| : ji I=l ~ N c , Ma — ~ 1/A C , we immediately obtain the result 



(3na ~ N*. (4.3) 

This surprising "superleading" behavior is the result of collective effects in the hedgehog 
wave function. In the orthodox approach to A c -counting, one should stop the analysis at 
this point, and conclude that as far as the magnetic polarizability is concerned, one is not 
close to the N c — > oo limit: according to Eq. (|4.3|) (3 should be much larger than a, which 
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contradicts the experiment. Using the physical values for various magnitudes in Eq. ( |4.2| ) 
we obtain the numerical value 



p NA ~ 12 x 1(T 4 fm 3 . (4.4) 

Note, that to compare fairly to experiment, a given model has to predict properly the 
value of f^NA- 

A related fundamental point is that the polarizabilities are defined in uj — > limits in the 
Compton amplitude, and in the physical world the resonances are at finite values of uj, well 
separated from the uj = region. In the large-iV c limit, the A-resonance occurs at uj = 0. 
Thus, the meaning of the polarizability is in fact lost, unless the uj — > limit is taken before 
the N c — > oo limit. In the linear response method, the limits are implicitly taken in the 
reversed order. 

In principle, a consistent analysis of the iV c -subleading physics is possible, but this would 
require a fully quantum-mechanical (not semiclassical) treatment, e.g. one could use the 
Kerman-Klein method [57]. Technically, this would be a tremendous effort, and of question- 



able merit in a simple effective theory such as hedgehog models. 

In order to be able to make some estimates, we relax the strict iV c -counting requirement 
in the rest of this section, and try to examine the subleading terms in N c in a more "flexible" 
approach. 



B. Sea-gull term 

In hedgehog models, the magnetic seagull term is related to the electric sea-gull term by 
the relation ( |2.17[ ). In enters at the level N*. Numerically, 

Psea-gull = ~U X 10" 4 fm\ (4.5) 

which would almost exactly cancel the (3na term, in apparent rough agreement with exper- 
iment. This is, however, not the full story at the level N c . In the next section we show, that 
the dispersive contributions to (3 arise of the same order as /3 sea -guih but have positive sign, 
and the desired cancellation is largely suppressed. 



C. Dispersive terms 

The calculation of the magnetic dispersive terms goes along the same lines as the calcu- 
lation of the electric polarizability. There are, however, a few important differences. Firstly, 
for the evei-TZ B 1 perturbation the pionic term in Eq. (|2.12|) leads to pionic sources in the 
linear response equations. Since these are long-ranged (proportional to 717J, we get strong 
dispersive effects. The sources are listed explicitly in Table |. For the K = ++ and 2 ++ 
cases we obtain equations of the form (II. 3. 8). The only difference is that in the K = ++ 
the valence quark eigenvalue, e, changes, and the resulting quark equation in (II. 3. 8) has 
the form 

{h - e)5q + - 5sq h - 2gMq h Z = j+, (4.6) 
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where Se is the shift in the quark eigenvalue. This shift may be viewed as a Lagrange 
multiplier ensuring the orthogonality of the shifted quark spinor, Sq + , and the hedgehog 
quark spinor, qv The numerical methods of treating the K = ++ and 2 ++ cases are 
straightforward. 

In the K = 1 ++ case, however, the perturbation excites the rotational zero-mode. This 
is a similar mechanism as in the case of the excitation of the translational zero mode in 
the E case - are the quantum numbers of the rotational mode in spin, or 

isospin (these are equivalent, since the hedgehog soliton is a K vn = ++ object). Thus, the 
hedgehog baryon placed in a constant magnetic field starts to rotate, and continues to spin 
up indefinitely. Quantum-mechanically, if time-dependent perturbation theory were used, it 
would correspond to transitions to higher spin and isospin states. 

There are a few subtle issues here, concerning the order of limits. Suppose we perform 
(exact) projection first (e.g. suppose we have done a calculation using the Kerman-Klein 
method |57|]), and have states with good angular momentum and isospin on which we apply 
the magnetic perturbation. Then, the rotational zero modes no longer appear, since they 
arise only if the solution breaks the symmetry of the lagrangian (Sec. Ill B in (II)). Thus, 
the appearance of the rotational zero mode in the K vn = 1 ++ is the effect of starting from 
the unprojected hedgehog. In fact, this problem is another manifestation of the noncom- 
mutativity of the u — > and the N c — > oo limits, described in Sec. (|1V A| ). One may hope, 
that in an improved treatment, the zero-mode physics is described by the N — A term, Eq. 
( |4.2| ). With this in mind, in order to estimate the size of the dispersive effects, we project 
out the divergent zero- mode part from the K = 1 contribution, and retain the physical part. 
The total result for (3 consists of /3na, Psea-guii, and fidisp. with the zero- mode contribution 
projected out. As mentioned before, this is not a consistent procedure, but it allows us to 
estimate the importance of various effects in the magnetic polarizability. 

The numerical values obtained for the K — 0, 1 and 2 parts of the B 1 magnetic pertur- 
bation are 3.6, —0.6 and 3.3 x 10~ 4 /m 3 , respectively. The total dispersive contribution to 
Pn is 

(3 disp = 6.5 x KT 4 /™ 3 , (4.7) 

The quark contribution is negative, and carries —15% of the total. 

The isoscalar magnetic perturbation, B°, which is a 1/N C effect, and comes entirely from 
the quarks, contributes 0.2 x 10 _4 /m 3 , which is negligible, as expected from N c counting. 
We also expect the neutron-proton splitting of magnetic polarizabilities, 5/3, to be small. It 
comes entirely from the quarks, since there is no magnetic analog of the electric term Q2.10| ), 
and is expected to be much smaller than 5a, which had a pionic contribution. 

At this point it might seem that adding up contributions ( j4.5| , [4~7j , fO| ) we still get a 
partial cancellation. However, the sea-gull and dispersive pieces are reduced if the pion field 
has the correct strength by a factor of ~ 2, as in the electric sea-gull case. The effects 
discussed in Sec. [V] further reduce these values, by another factor of 2. These effects cause 
the cancellation to disappear, and we are left with a large value of (3n- As stressed before, 
rigorous estimates could only be done in a calculation beyond linear response. 
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V. ROLE OF THE A IN HEDGEHOG MODELS AND CHIRAL PERTURBATION 

THEORY 



In this section we compare the hedgehog model predictions to the results obtained in 
chiral perturbation theory xPT . At first, it may seem awkward, since the two methods are 
based on two different limits: N c — > oo, and m w — > 0, which are known not to commute, 



and give different results for various observables ||20|| . However, as we show in Ref. ||58|| , for 
quantities which are divergent in the chiral limit as m' 1 , there is a simple connection between 
the hedgehog predictions and the X-PT predictions for scalar-isoscalar and vector-isovector 
quantities (electromagnetic polarizabilities ot^ and (3n are scalar-isoscalar). 

Let us evaluate the polarizabilities using our methods in the chiral limit. For ajv, the 
dominant part comes from the sea-gull term. In the chiral limit the quantity diverges as m^ 1 , 



and all of the contribution comes from the pionic tail, which has the form ( |3.2| ). Evaluating 
integral ( J2.17[ ) with the profile (^] we obtain 

aN = 32^MW (5 - 1} 

Similarly, including both the sea-gull and dispersive pieces for the magnetic polarizability, 
we obtain in the chiral limit 

= r /m2 • ( 5 - 2 ) 

These expressions are exactly a factor of 3 larger than the analogous expressions following 
from xPT fl5"5" |. This difference comes from the different treatment of the A isobar. In xPT 



the A contributions in pionic loops (Fig. ||) are not included, since it is implicitly assumed 
that the N — A mass splitting is much larger than the pion mass, and consequently loops 
with the A do not contribute to the leading singularity. They are counted as effects of order 
log(m n ), 1, .... In hedgehog models, on the contrary, the N — A mass splitting is a N~ l 
effect, much smaller than the pion mass, which is of the order N®. Therefore, hedgehog 
models include the A on equal footing with the nucleon. As shown below, spin and isospin 
Clebsch factors account for the difference by a factor of 3 between Eqs. ( |5.1| , |5.2|) , and xPT 
predictions. 

One may ask, how the physics of pionic loops, such as in Fig. |], is present in hedgehog 
models. After all, the treatment of the pion field in hedgehog models is classical. To 
demonstrate how hedgehog expressions in the chiral limit may be viewed as hadronic loop 
diagrams such as in Fig. f|, let us evaluate this diagram in the chiral limit. First, we perform 
a calculation of the self-energy due to an insertion of the electric sea-gull interaction in the 
pion loop. In momentum space, we obtain 

9A " 2 f f A r al5 (/k + y q )S F {p-k) 



V 2F/ J (2tt) 4 7 (2itY 
rM/k ~ \/q)D(k - \q)V{q)D{k + \q)T ah , (5.3) 

where Sp is the Feynman propagator of the nucleon, D(k) = i / '(k 2 —m^+ie) is the propagator 
of the pion field, V(q) = ^e 2 (27r) A 5 4: (q)^(E ■ Vq) 2 is the Fourier transform of the electric 
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sea-gull interaction in the coordinate space, V(r) = he 2 (E ■ r) 2 , and T ab = 5 a b — S a3 5b3- Now 



we proceed as follows: we first carry out the integral over k° in (|5.3| ). For the leading chiral 
singularity piece, m~ l , the contribution comes from the poles in the pionic propagators. The 
poles in the nucleon propagators contribute to less singular terms. We can then perform the 
nonrelativistic reduction of the nucleon propagator and the pion-nucleon vertices, and take 
the expectation value of S* in positive energy spinors, in order to extract the energy shift. 
Transforming back to coordinate space, we obtain the expression 

a = e 2 J d 3 xcf> a a Svm ir)(f> a b syrn ir)T ab , (5.4) 

which has the form of the expression obtained in hedgehog models. 

We now go back to the question of the factor of 3. When collective coordinates are 
introduced, expression ( |5.4|) becomes 

a = e 2 [d 3 x {<t> asym - x c) 2 . (5.5) 



In Sec. |111 A| , to obtain the sea-gull contribution to the oln we evaluated the collective matrix 
of expression ( |5.5|) in nucleon collective wave functions. Now we repeat the procedure, but 
we write the resulting matrix element of the integrand in Eq. (|5.5|) as 



]T£(iV|(<r^ x c) a \i)(i\(cf> as ^ x c) a \N), (5.6) 

i a 

where \i) is a collective nucleon or A intermediate state. If the sum over i is unrestricted, 
then we just recover our previous expressions, (|5.1| , |5.2|) . If, however, we restrict i to run 
only over nucleon collective states, then using Eq. (BIO) in (II), we obtain a result smaller 
by a factor of 3, and this result agrees with the xPT prediction. Conversely, had the xPT 
calculation included the A in diagrams of Fig. (|, with Ma = Mn, it would predict a result 



3 times larger than quoted in Ref. [[55 
Let us introduce 

d={M A - M N )/m n . (5.7) 

These two cases, hedgehog models and xPT ^ correspond to two limits: d — > 0, and d — > oo. 
In nature, d ~ 2, which is between the two limits, and we do not have the separation of the 
pion mass and Ma — M N scales. In this case it seems most appropriate to treat both scales 
as small, and keep d as an unconstrained parameter. This is the spirit of the approach of 
Refs. To estimate the contribution of the A at the physical value of d, we evaluate 



the diagram of Fig. |], starting from expression ( |5.3| ), with Sf describing the A propagator 
with the physical mass, and with the vertices modified appropriately. The m" 1 contribution 
is easily obtained, since, as in the nucleon calculation presented above, we can perform 
the non-relativistic reduction. The result for the ratio of the A to nucleon contribution in 
diagram |4] is 

^ = 25(d), (5.8) 

where the 2 comes from Clebsch factors, and the "mass suppression" function has the simple 
form 
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, , 4 f Arctan (J±=£) for d < 1 

5(d) = -< V Y 1 £H< , - 5.9 

W 7T | Arctanh / v 7 ^ 77 ! for d > 1 V ; 

For degenerate N and A, S(d = 0) = 1. In the limit of d — > oo, S(d) ~ logd/d. For 
the physical value of d we find 5(2.1) = 0.47, and from Eq. (|5.8p we find, that the A 
contribution in pionic loops of Fig. [| to the electric polarizability is roughly equal to the 
nucleon contribution. The results for the magnetic polarizability are analogous. 

This important feature of large A contributions in pionic loops (for scalar-isoscalar op- 
erators) is also present in other quantities []5"8|1 . In the above discussion we have tacitly 
assumed that the ratio g w NA/ ' 9-kNN = 3/2 p0| , |5l| , as predicted by hedgehog models. Exper- 



imental numbers agree with this prediction to within a few percent. The correction may be 
introduced to expression (|5.8|) |58[ . 



We note that the above analysis and the comparison of the hedgehog model predictions 
and the xPT predictions can be done only for observables which do not depend on cranking 
(independent of the cranking frequency, A). Observables, which do depend on cranking do 
not have the same chiral singularities, e.g. the electric mean squared radius diverges as m^ 1 
in hedgehog models, and as logm 71 " in ;^PT. The same is true of the splitting of electric 
polarizabilities, 5a. In this case the issue of nonocommutativity of the large- N c and chiral 
limits is much more complicated. 

The results of this section indicate how one should try to improve the hedgehog pre- 
dictions by subtracting the amounts by which these models overestimate the leading chiral 
singularity terms. There is uncertainty in such a non-rigorous procedure, but in our view it 
is required by physics. As already mentioned at the end of Sec. |IV A| , the correct treatment 
of the A could be done in a quantum calculation with the inclusion of the iV c -subleading 
terms. 

On the other hand, the predictions of X-PT will be largely effected by the diagrams of 
Fig. |] with A intermediate states. For polarizabilities, at the leading singularity level, these 
diagrams are as important as the diagrams with the nucleon intermediate states. 



VI. EFFECTS OF THE PIONIC SUBSTRUCTURE 

The minimal substitution prescription in an effective lagrangian cannot produce all in- 
teractions of a hadronic system with the electromagnetic field. It is possible to write down 
terms which are by themselves gauge-invariant, and thus not obtainable by gauging la- 
grangian Q2.1p. In this section we focus on terms £ 9 and £ i0 of reference jSTfl: 

C 9 = -iL 9 Tr[^ v D^UD u U^ + T R V D^ D V U\, (6.1) 
C 10 = L 10 Tr[^ u U^ R U% (6.2) 

where Tfy R are the left and right chiral field strength tensors. For the case of electromagnetic 
field we have Fb, R = e^r^F^ = e^r 3 (d^A u — d v A^). In the linear a-model, U corresponds 
describes the chiral field, U = F^{a + it • 7r). The constants L 9 and L w can be expressed 
through measurable quantities, namely the pion electric mean square radius, (r 2 )^, and the 
pion polarizability a 71 = —(3 : 
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F 2 (r 2 V 

U = ^£ (6.3) 

L 10 = -L 9 = m ^ a \ (6.4) 

In the chiral perturbation theory treatment, one considers pionic loops, and through renor- 
malization L g and L w acquire chiral logarithms. In our approach, we simply treat the terms 
( |6.2| ) in the mean-field approximation, replacing the meson field operators by classical fields. 
As before, the E and B fields are constant, and we find from the £9 term 

d 3 xC q = 4e 2 L q F~ 2 f d 3 xF iu AJc x Tz h ) ■ (c x 8^,) 



= 2L 9 (E 2 - B 2 ) J d 3 x(c x 7z h ) 2 , (6.5) 

where i runs over spatial values. The £10 term leads to 

J d 3 xC 10 = 4e 2 L 10 (E 2 - B 2 )F~ 2 

J d 3 x ((cxtt^-K^ + t^)). (6.6) 

The last term under the integrand is canceled by the vacuum subtraction, which is implicit. 
It becomes — {a\ + 7r|) + a 2 ac = —(u 2 + 7r^) + F 2 , which is zero in the nonlinear sigma 
model, but also vanishingly small in our case due to proximity of the hedgehog solution to 
the chiral circle. We can thus drop this term in our estimate. Using (|6.4| , II. 3. 21) we read 



off from ( |6.5| , 3.6 ) the following contribution to the nucleon polarizabilities: 

a N = — p N = m T U a . (6.7) 

The physical interpretation of this contribution is clear. The pion, having electromagnetic 
structure, is polarizable. Since the nucleon is surrounded by a pion cloud, this pion polar- 
izability results in additional polarizability of the nucleon. Note the opposite signs of the 
electric and magnetic polarizabilities in ( p.7|) , reflecting the fact that ~6F = — /3 7r . Also notice, 
that since oF is of order 1 in iV c -counting, the contributions |6.?1 are of order iV c . The terms 
( |6.2| ) do not lead to additional neutron-proton splitting of polarizabilities 

We can call N n = m T 9 mes in |6.7| the "number of pions" in the nucleon seen in the 
Compton scattering process. Numerically, N n ~ 0.5 in our model. Using the relation of the 
moment of inertia to the N — A mass splitting, we can write 

where the last factor is the fraction of the total moment of inertia carried by the pion. This 
quantity is ~ 1 in hedgehog models: 60% in our model with quarks, 100% in the Skyrmion. 
We thus have a quasi-model-independent result N n ~ 0.5 — 0.7. 

The value of a n can be determined experimentally, however existing experimental data 
)1| do not seem reliable, and are in contradiction with a low-energy theorem due to Holstein 



49fl , which gives a n = 2.8 x 10 4 /m 3 . With this value we get 
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ajj = -(31 ~ 1.3 x 10 _4 /m d , (6.9) 

which is a few times smaller compared to the minimal substitution terms, but non- negligible, 
especially for the magnetic case, where we expect cancellations to occur. If experimental 



numbers for a n were used |6l] , then a three times larger result would follow. 

Other non-minimal substitution terms can also be considered, but there is no knowledge 
of the phenomenological low-energy constants which have to be introduced, and predictive 
power is lost. 



VII. OTHER MODELS 



The analysis of this paper can be straightforwardly applied to other chiral models. In 
purely mesonic Skyrmions ||21|| , the RPA approximation corresponds to linearizing the small 



fluctuation equations, hence the non-linearity constraint is not imposed at the quantum level. 
Appropriate linear response equations can be derived along the lines of App. B in (II). The 
role of the isoscalar source is carried by the topological (Goldstone-Wilczek) current, and 
effects of the higher order terms in the equations are present. 

In purely quark models (NJL) ||32|| , upon minimal substitution, the lagrangian in presence 
of the electromagnetic field A^, and after introducing collective degrees of freedom, becomes 

Snjl(A m ) = -iTrlog (i/d - g(a + ry 5 T • tt) - ±7 T • A 

+ \/A(N~ 1 + r • cj) - (vac), (7.1) 



where (vac) denotes the vacuum subtraction, and an NJL cut-off is understood. The trace 
accounts for the occupied (valence) levels. This expression may be expanded to second-order 
in A^, and expressions for polarizabilities can be easily derived. The dispersive pieces, as 
in the model treated in this paper, lead in case of even- TZ interactions to distortions in the 
profile functions a and 7r. This is likely to result in numerical complications, in particular in 
cases with zero-mode excitations. One can use gradient expansion techniques instead of 
solving the model exactly. Then, for example, the isovector electric responce of the quarks 
from the Dirac sea generates the pionic sea-gull contribution. 

All our generic hedgehog model conclusions described in this paper hold for these, and 
other hedgehog models. 



VIII. CONCLUSION 

The main results obtained in this paper can be summarized as follows: 

1) We show the A^-counting rules for electromagnetic polarizabilities in hedgehog models. 
For electric polarizabilities the basic experimental pattern is reproduced (a^ ~ N c > 5a ~ 
1/N C ). For the magnetic case the rules show the inapplicability of linear response (fix ~ N%). 

2) Dispersive terms lead to deformation of hedgehog solitons, and the resulting con- 
tributions to polarizabilities enter at the same iV c -level as the sea-gull contributions. In 
the magnetic case, we expect that the presence of dispersive terms leads to large positive 
contributions to (3n- There is no cancellation mechanism which can bring (3 down to the 
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experimental value. In fact, all model calculations with the A degree of freedom have this 
problem, since the paramagnetic N-A term is large, and there is no simple mechanism to 
cancel this effect. 

3) Hedgehog models provide a mechanism of splitting of the neutron and proton polar- 
izabilities. An explicit calculation gives reasonable numbers, and the sign of 5a = a n — a p 
is expected to be positive in models with pionic clouds. 

4) Concerning numerical results in hedgehog models, because of the sensitivity of the 
results to the pion tail, the value of qa in a model should be well reproduced. Also, for the 
magnetic case, good prediction for [ina, as well as Ma — M N , is necessary to reproduce the 
N — A paramagnetic term. 

5) The A resonance plays an important role. Hedgehog models largely overestimate 
these contributions, since they neglect the effects of N — A mass splitting in projection. 
We show how to estimate the A effects in pionic loops in a modified chiral perturbation 
theory, and find these effects are at the level of 100 % in calculation of the electromagnetic 
polarizabilities. One should also try to improve the hedgehog predictions by subtracting the 
amounts by which these models overestimate the leading chiral singularity terms. 

6) The effects of non-minimal substitution terms Cg and C±o |5(J in the effective la- 
grangian enter at the level of 1 — 2 x 10~ 4 /m 3 . 
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FIGURES 



FIG. 1. Dispersive contributions to the Compton scattering amplitude. In our treatment, Q is 
the RPA (linear response) propagator. 

FIG. 2. Radial densities of the leading- iV c contributions to ajy (in units of 10~ 4 /m 2 ): Sea-gull 
(solid line), and quark E 1 contribution (dashed line). The dotted line shows the sea-gull term 
evaluated with the asymptotic profile, Eq. ( |3.2| ). 

FIG. 3. Radial density of the pion contribution to the splitting of the neutron and proton 
electric polarizabilities, 5a (in units of 10 _4 /m 2 ). 

FIG. 4. Hadronic one-pion-loop diagrams giving the leading chiral contribution to nucleon 
polarizabilities. The vertex corresponds to the sea-gull interaction e 2 / d?x (E ■ x)(c X 7r) 
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TABLES 



TABLE I. K PR classification of various dispersive perturbations, and sources in the corre- 
sponding linear response equations (App. A in (II)). Quantum numbers (L,A) label quark shifts, 
L n and L a are the orbital numbers of meson shifts. Null entry in columns for L n and L a means 
that the fluctuation does not arise. For cases with zero modes (electric isoscalar and magnetic 
isovector, K=l), the sources listed are twice the jx, jy or jz sources from equations in App. A in 
(II). For other cases they are the sources entering Eqs. (11.3.7,11.3.8). Radial functions Gh and Fh 
are the upper and lower components of the hedgehog valece spinor, and in the hedgehog pion 
profile. See (II) for details. 
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TABLE II. Electric polarizability. 


The model predictions for a at 


are expected to be largly 


reduced by effects disussed in the text. 
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